arXiV: astro-ph/0211412v2 17 Feb 2003 


Draft version February 2, 2008 

Preprint typeset using style emulateapj v. 11/12/01 


HIGH-RESOLUTION CHANDRA SPECTROSCOPY OF TAU SCORPH: A NARROW-LINE X-RAY 

SPECTRUM FROM A HOT STAR 

David H. Cohen and Genevieve E. de Messieres 
Department of Physics and Astronomy, Swarthmore College, Swarthmore, PA 19081 
cohen@astro.swarthmore.edu, gdemessl@swarthmore.edu 

Joseph J. MacFarlane 

Prism Computational Sciences, 455 Science Dr., Madison, WI 53711 
jjm@prism-cs.com 

Nathan A. Miller 

Department of Physics and Astronomy, 105 Garfield Ave., University of Wisconsin, Eau Claire, Eau Claire, 

WI 54702 
millerna@uwec.edu 

Joseph P. Cassinelli 

Department of Astronomy, 475 N. Charter St., University of Wisconsin, Madison, Madison, WI 53706 

cassinelli@astro.wisc.edu 

Stanley P. Owocki 

Bartol Research Institute, University of Delaware, Newark, DE 19716 
owocki@bartol.udel.edu 

AND 

Duane A. Liedahl 

Lawrence Livermore National Laboratory, Livermore, CA 94550 
liedahll@llnl.gov 
Draft version February 2, 2008 

ABSTRACT 

Long known to be an unusual early-type star by virtue of its hard and strong X-ray emission, t Scorpii 
poses a severe challenge to the standard picture of 0-star wind-shock X-ray emission. The Chandra 
HETGS spectrum now provides significant direct evidence that this BO.2 star does not fit this standard 
wind-shock framework. The many emission lines detected with the Chandra gratings are significantly 
narrower than what would be expected from a star with the known wind properties of t Sco, although 
they are broader than the corresponding lines seen in late-type coronal sources. While line ratios are 
consistent with the hot plasma on this star being within a few stellar radii of the photosphere, from at 
least one He-like complex there is evidence that the X-ray emitting plasma is located more than a stellar 
radius above the photosphere. The Chandra spectrum of t Sco is harder and more variable than those 
of other hot stars, with the exception of the young magnetized O star 6^ Ori C. We discuss these new 
results in the context of wind, coronal, and hybrid wind-magnetic models of hot-star X-ray emission. 

Subject headings: circumstellar matter - stars: early-type - stars: mass loss - X-rays: stars 


1. INTRODUCTION 

The surprising discovery of X-ray emission from OB 
stars by the Einstein satellite in the late 1970s (Hamden 
et al. 1979; and anticipated by UV superionization ob¬ 
servations [Cassinelli & Olson 1979]) is now a 25-year-old 
mystery. Because magnetic dynamos and the associated 
coronal X-rays were not believed to exist in hot stars, al¬ 
ternative models have been sought since the discovery of 
X-rays from hot stars. Most of these center on some type 
of wind-shock mechanism that taps the copious energy in 
the massive radiation-driven winds of hot stars (Lucy & 
Solomon 1970; Castor, Abbott, & Klein 1975, hereafter 
CAK). Wind-shock models, especially the line-force in¬ 
stability shock model (Owocki, Castor, & Rybicki 1988; 
Feldmeier, Puls, & Pauldrach 1997; Feldmeier et al. 1997), 
have become the favored explanations for the observed OB 
star X-ray properties. 

However, various pieces of evidence have recently 
emerged indicating that magnetic fields exist on some hot 


stars (Henrichs et al. 2000; Donati et al. 2001, 2002). 
There is also mounting evidence that the standard wind- 
shock scenario is inadequate for explaining the observed 
levels of X-ray production in many B stars (Berghofer & 
Schmitt 1994; Cohen, Cassinelli, & MacFarlane 1997a; Co¬ 
hen, Cassinelli, & Waldron 1997b), and perhaps in some 
O stars as well (Schulz et al. 2001; Waldron & Cassinelli 
2001; Miller et al. 2002), leading to a revival in the idea 
that some type of dynamo mechanism might be operat¬ 
ing on hot stars. Theoretical investigations of magnetized 
OB stars have followed (Charbonneau & MacGregor 2001; 
MacGregor & Cassinelli 2002), including studies of hy¬ 
brid wind-magnetic models (Babel & Montmerle 1997; ud- 
Doula & Owocki 2002; Cassinelli et al. 2002). 

Discriminating between the two broad pictures of X-ray 
emission in hot stars (magnetic coronal vs. wind shock) 
has been very difficult during the era of low-resolution X- 
ray spectroscopy. But with the launch of Chandra in 1999, 
a new, high-resolution window has opened, and we can now 
apply quantitative diagnostics in our attempts to under- 
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stand the high-energy mechanisms operating in the most 
massive and luminous stars in our galaxy. 

Coronal heating mechanisms are complex and, even in 
our own Sun, not fully understood. Their possible appli¬ 
cation to hot stars is certainly not clear, even in those 
hot stars for which magnetic fields have been detected. 
However, we can assume that any type of coronal X-ray- 
production mechanism involves the confinement of rela¬ 
tively stationary hot plasma near the photosphere. This 
is in contrast to wind shocks, which, whatever their na¬ 
ture (and many different specific models have been pro¬ 
posed), involve low-density hot plasma farther from the 
photosphere and moving at an appreciable velocity. This 
broad dichotomy naturally leads to several possible ob¬ 
servational discriminants that might be implemented us¬ 
ing high-resolution X-ray spectroscopy. These include the 
analysis of line widths (or even profiles), which could re¬ 
flect Doppler broadening from the bulk outflow of winds 
(which in normal OB stars have terminal velocities of up 
to 3000 km s“^) and emission-line ratios that are sensitive 
to density, the local radiation field, or (unfortunately) in 
some cases, both. There are numerous other spectral diag¬ 
nostics available, including temperature-sensitive line ra¬ 
tios, spectral signatures of wind absorption, time variabil¬ 
ity of individual lines, and differential emission measure 
(DEM) fitting (multiline plus continuum global spectral 
modeling of X-ray emission from plasma with a continuous 
distribution of temperatures). All of these can be useful 
for determining the physical properties of the hot plasma 
on OB stars, and, ultimately for constraining models, but 
none can provide the same level of direct discrimination 
between the wind-shock and coronal scenarios as line pro¬ 
files and ratios. 

In this paper we report on the very high signal-to-noise, 
high-resolution Chandra HETGS spectrum of the BO.2 V 
star r Sco. We chose to observe this particular hot star 
because it is one of the most unusual single hot stars, 
with indications already from various low-resolution X- 
ray observations that its X-ray properties are extreme and 
present a severe test to the standard wind-shock scenario. 
As such, it makes an interesting contrast to the five single 
O stars that have been observed at high resolution with 
Chandra and XMM up to this point. In fact, some of 
the unusual properties of r Sco (discussed below in §3) 
have already inspired the development of alternate models 
of X-ray production, notably including the “clump infall” 
model of Howk et al. (2000). In this scenario overdense re¬ 
gions condense out of the wind, and the reduced line force 
due to the increased optical depth in the clumps results 
in these clumps falling back toward the star, generating 
strong shocks via interactions with the outflowing mean 
wind. 

2. OBSERVATIONS 

The data presented in this paper were taken between 
2000 September 17 and 18, over two different Chandra or¬ 
bits, with an effective total exposure time of 72 ks. We 
used the ACIS-S/HETGS configuration and report on the 
medium-energy grating (MEG) and high-energy grating 
(HEG) spectra, although we concentrate primarily on the 
MEG 1st order spectrum, which contains a total of 19348 
counts (while the HEG 1st order spectrum has 5230 to¬ 


tal counts). We extracted and analyzed the spectra with 
CIAO tools (Version 2.1.2), including the SHERPA model¬ 
fitting engine. We used the GALDB libraries. Version 2.6, 
and the OSIP files from 2000 August to build auxiliary 
response files and grating response matrices. We fit line- 
profile models using SHERPA, but astrophysical modeling 
of the spectral data took place outside of the CIAO envi¬ 
ronment with a custom non-LTE spectral synthesis code, 
described in section 5.2. 

3. THE STAR, T SCO 

A member of the upper Sco-Cen association and a Mor¬ 
gan and Keenan system standard, the BO.2 V star r Sco 
has been the subject of much study over the past several 
decades. At a distance of 132 ± 14 pc (Perryman et al. 
1997), this niy = 2.8 star is closer to the Earth than any 
hotter star. It has an effective temperature of 31,400±300 
K and log g = 4.24 ± .03 (Kilian 1992). UV observa¬ 
tions show evidence for a wind with a terminal velocity 
of roughly 2000 km s“^ (Tamers & Rogerson 1978). The 
mass-loss rate is difficult to determine because of the un¬ 
certain ionization correction, but it seems consistent with 
the theoretically expected value of M « 10“® Mq yr 
(calculations based on Abbott 1982). Published values 
based on UV absorption lines (e.g.. Tamers & Rogerson 
1978) and a recently determined upper limit based on IR 
emission lines (Zaal et al. 1999) are roughly half this value. 
These wind properties contrast with those of the prototyp¬ 
ical early O supergiant C Pup, which has a mass-loss rate 
at least 100 times larger (because of its significantly larger 
luminosity). Chandra HETGS observations of C Pup were 
recently reported by Gassinelli et al. (2001), and clear 
signatures of the stellar wind (broadening and absorption 
effects) are seen in those data. 

There are several unusual properties of r Sco that might 
be relevant for our attempts to understand its X-ray prop¬ 
erties: (1) It has a quite highly ionized wind, since the 
star is the coolest to show O VI in its COPERNICUS 
spectrum (Tamers & Rogerson 1978). (2) The wind is also 
stronger than normal for the star’s spectral type (Walborn, 
Parker, & Nichols 1995). (3) It has a very low Vrots'mi of 
10 ± 5 km s“^ (Wolff, Edwards, & Preston 1982). (4) It 
shows unusual redshifted absorption of about 200 km s“^ 
in some of the higher ionization UV lines (which also show 
much more extended blue absorption, as one would expect; 
Tamers & Rogerson 1978). (5) It is extremely young; per¬ 
haps being only 1 million years on the main sequence (Kil¬ 
ian 1994). (6) It shows unusually strong photospheric tur¬ 
bulence, as seen in optical and UV absorption lines (Smith 
& Karp 1979) and in the IR (Zaal et al. 1999). 

In addition to these properties, it is already known from 
Einstein, ROSAT, and ASCA measurements that r Sco 
has an unusually hard X-ray spectrum and a high X-ray 
luminosity (Gohen et al. 1997a, 1997b). It is this fact, 
along with the redshifted UV absorption (item 4 above), 
that prompted Howk et al. (2000) to develop their model 
of infalling clumps for r Sco. 

Many of these properties have been cited as evidence 
that r Sco may not fit into the usual wind-shock picture 
of hot star X-ray production, or at least that its wind 
emission mechanism may be modified by the presence of 
a magnetic field or perhaps involve some other, less stan- 
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dard, wind-shock model. We address these issues in detail 
in the next sections. 

4. THE DATA 

The overall quality and richness of the Chandra HETGS 
data can be seen in the MEG hrst-order spectrum shown 
in Figure 1. There are numerous lines present in this high 
signal-to-noise spectrum, and they are much better sepa¬ 
rated than in any of the other hot stars thus far observed 
(because of the small wind broadening seen in r Sco com¬ 
pared to these stars). In Figure 2 we show a portion of the 
MEG spectrum, centered on the neon Lya line, compared 
with data from C Pup and Capella, a typical O supergiant 
wind source and a late-type coronal source, respectively. 
It is immediately obvious that t Sco is superhcially more 
similar to Capella than to C Pup. 

To quantify the properties of the spectrum, we fitted 
every strong line with a simple model (Gaussian plus a 
polynomial for the nearby continuum), providing values of 
the line flux, characteristic width, and central wavelength. 
We employed the appropriate spectral response matrices 
in these fits, so the derived line widths are intrinsic and 
do not include instrumental broadening. When there was 
sufficient signal, we fit the MEG and HEG data simultane¬ 
ously. These results are summarized in Table I, with the 
several strongly blended spectral features listed in Table 
2. The astrophysical analysis of these data is presented in 

. . . . 

Time variability properties of r Sco are also of interest in 
discriminating between the two primary X-ray-production 
mechanisms. In Figure 3 we show a lightcurve formed of 
all the photons in the -l,-2,-3,-|-3,-|-2,-|-l orders of both the 
MEG and the HEG. The bins are 1000 seconds in length. 
The data are not consistent with a constant source (at 
P > 99.99 %), but no large flare-type events, with a sud¬ 
den increase in flux followed by a gradual decline back to 
the basal level, are seen. The standard deviation of the 
data is slightly less than 10 % of the mean. Only one indi¬ 
vidual line feature, the oxygen Lya line near 19 A, shows 
evidence of significant variability. Its lightcurve is shown 
in Figure 4. We tested other strong lines, and only the 
Fe XVII line at 15.01 A shows even marginal variability. 

5. THE X-RAY DIAGNOSTICS AND ANALYSIS 
5.1. Line Widths 

The emission lines of r Sco are obviously narrow com¬ 
pared to the same lines observed in hot-star wind X-ray 
sources like C Pup, but how intrinsically narrow are they? 
Might they be no broader than thermal broadening would 
dictate? To answer these questions, we fitted Gaussian line 
models, convolved with the instrumental response func¬ 
tion, so the width of the model Gaussian reflects the actual 
width of the emitted line. We do not know what actual 
line prohle shape is warranted, but because the Gaussian 
model is simple and relatively general and the observed 
line broadening is not dramatic, we did not attempt to fit 
other profile shapes. 

The line widths-half-width at half maximum (HWHM)- 
listed in Table 1 generally exceed by a factor of about 3 
those expected from thermal broadening (assuming a tem¬ 
perature equal to that at which the emissivity for a given 
line peaks). In addition, fits we performed to Gapella and 


AB Dor (both late-type coronal sources) using the same 
procedures found narrower lines (consistent with thermal 
broadening) than we found for t Sco. In Figure 5 we show 
the Ne X line of r Sco in the MEG, along with an intrinsi¬ 
cally narrow line (Gaussian model with ct = 0; i.e. a delta 
function) convolved with the instrumental response. The 
data clearly are broader than the model. We also note 
that some of the features in the spectrum are doublets. 
We generally htted these with single-Gaussian models. In 
each case in which we fitted a doublet with a single-line 
model, the separation of the two components of the dou¬ 
blet was far less than the typical (and derived) line width. 
Thus the derived line widths using this procedure are only 
marginally elevated because of the blended doublet. 

The velocities of most of the lines in the HETGS spec¬ 
trum of r Sco are a few hundred km s“^ (HWHM), which, 
while greater than the thermal velocity, is signihcantly less 
than the wind terminal velocity. The X-ray line velocities 
seen in t Sco are thus larger than those in coronal sources 
and smaller than those in wind sources, a trend that is 
shown in Figure 6, and importantly, the velocities derived 
for the wind of t Sco are significantly less than what would 
be seen if the X-ray-emitting gas had the same velocity 
distribution as the UV-absorbing gas. 

The mean HWHMs of the strongest, unblended lines are 
roughly 300 km s“^, or 15% of the terminal velocity of the 
wind of T Sco (see Fig. 7). There is a slight trend of de¬ 
creasing line width with wavelength, which is the opposite 
of what is seen in the O stars (Kahn et al. 2001; Waldron & 
Gassinelli 2001; Gassinelli et al. 2001; Miller et al. 2002). 
In the O stars, the trend of increasing line width with 
wavelength can be understood in terms either of temper¬ 
ature stratification (lower temperature plasma farther out 
in the wind, traveling at higher velocities) or an absorption 
effect (low energies are preferentially attenuated, so these 
lines tend to be broader because only photons from the 
outer wind, where velocities tend to be high, can escape). 
The lack of such a trend in r Sco might be evidence of 
the lack of any wind continuum absorption effects on the 
observed spectral lines. The wind is expected to be opti¬ 
cally thin down to r < 1.5 i?* near 0.5 keV and effectively 
completely optically thin at photon energies of 1 keV and 
above. 


5.2. Radiation-Sensitive Line Ratios 

The He-like forbidden-to-intercombination line strength 
ratio is known to be sensitive to the local mean inten¬ 
sity of the UV radiation field (Gabriel & Jordan 1969; 
Blumenthal, Drake, & Tucker 1972; Waldron & Gassinelli 
2001; Kahn et al. 2001). If the photospheric flux is strong 
enough, then electrons can be radiatively excited out of 
the metastable upper level of the forbidden line (^^i) to 
the upper level of the intercombination line (^^ 1 , 2 )) weak¬ 
ening the forbidden line and strengthening the intercombi¬ 
nation line. Thus a measurement of the f/i ratio [using the 
more standard spectroscopic notation, z/(x y)], poten¬ 
tially yields the radius of the X-ray emitting plasma (via 
sensitivity of the mean intensity to the dilution factor). 
This assumes that we know or can ignore limb darkening 
and other geometrical effects and that we know the UV 
luminosity of the star at the relevant wavelengths. 
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More traditionally, this line ratio has been used as a 
density diagnostic in coronal sources (where, because of 
the low temperature of the stars, very little UV radiation 
is present). Of course, in hot stars, this line ratio is still 
sensitive to density (via collisional excitation between the 
same two levels described above). However, the densities 
that would be required to explain the small observed f/i 
values are very large, while the radiation field requirements 
are much more reasonable. 

In order to extract information from these line ratios, 
we model the oxygen, neon, magnesium, and silicon ex¬ 
citation/ionization using NLTERT, a non-LTE radiation 
transport and statistical equilibrium code (described in 
MacEarlane et al. 1993; MacFarlane, Cohen, and Wang 
1994). We use a Kurucz model atmosphere, constrained 
by observational data in the spacecraft UV (relevant for 
O and Ne), but unconstrained by direct observation in 
the far-UV (Mg) and extreme-UV (EUV) (Si and S). We 
note that the model atmosphere fluxes are sensitive to the 
adopted temperature only in the EUV, in which a change 
of 1000 K can lead to about a factor of 2 change in emer¬ 
gent flux. 

The Hea line complexes are shown in Figure 8. When 
the metastable level is neither collisionally nor radiatively 
destroyed (i.e., the “low-density limit”), the forbidden 
lines tend to be about a factor of 3 stronger than the in¬ 
tercombination lines. We can see in the figure that the 
forbidden line is heavily affected in oxygen and neon, less 
so in magnesium, and less still in silicon. The S XV data 
are marginal, but are shown for the sake of completeness. 
In Figure 9 we show the results of our model calculation 
for He-like magnesium. Note that we model the effects of 
both radiation and collisions (i.e., showing both the ra¬ 
dial sensitivity and the density sensitivity) in these cal¬ 
culations. Similar calculations were carried out for the 
He-like features for the other four species. The results of 
these calculations and the constraints they imply for the 
location of the X-ray emitting plasma are summarized in 
Table 3. The quoted errors on the determinations of the 
radii of line formation include statistical errors only. The 
uncertainties in the atmosphere model and the associated 
UV luminosity will certainly affect the results for Si XIII 
(and S XV) more than for the other three line ratios. For 
example, if the flux were actually twice as high as in the 
model atmosphere, then the lower limit to the radius of 
line formation would move out from r = 1.1 to almost 
1.5 R*. 

We also point out that, where other lines were present 
(primarily iron lines near the neon complex), we fitted 
them at the same time we fitted the He-like features, and 
upper limits to the nondetected forbidden lines were de¬ 
termined by fitting a Gaussian model with its centroid 
position fixed at the laboratory wavelength and its width 
fixed at the value derived for the corresponding He-like 
resonance line. 

In general, it can be seen that most of these line ratio 
diagnostics are consistent with the X-ray emitting plasma 
being within several stellar radii of the photosphere. The 
oxygen and sulfur complexes do not really provide inter¬ 
esting constraints, as they are consistent with the X-ray 
emitting plasma being either right above the photosphere 
or far out in the wind flow. The neon ratio implies that the 


plasma giving rise to this complex is within a stellar ra¬ 
dius of the photosphere, but confusion and blending with 
nearby iron lines may make this result somewhat less firm 
than the formal errors imply. The strongest constraints are 
provided by the Mg XI and Si XIII features. The Mg XI f/i 
ratio requires the X-ray emitting plasma be at least 1.6 i?* 
above the photosphere (but no more than 1.9 R/. The 
Si XIII f/i ratio indicates that the hotter plasma giving 
rise to this feature is closer to the star, a trend that is also 
seen in the O stars observed with Chandra (Waldron & 
Cassinelli 2001; Cassinelli et al. 2001; Miller et al. 2002). 

Finally, we note that the intercombination line fluxes 
exceed the resonance line fluxes in the case of oxygen and 
neon. While the intercombination line strength is gener¬ 
ally expected to be less than that of the resonance line 
in collisional plasmas, the observed G ratios {{i + f)/r) 
for these complexes are not inconsistent with a colli¬ 
sional plasma, since the enhanced intercombination line 
flux comes at the expense of the weakened forbidden line 
(Pradhan & Shull 1981). 

5.3. Density-Sensitive Line Ratios 

While the helium-like forbidden-to-intercombination ra¬ 
tios are sensitive to both density (via the collisional de¬ 
struction of the level) and the local radiation field 
(through photoexcitation from the same level), we have 
argued in §5.2 that in hot stars like r Sco, the ratio is 
dominated by radiative effects. This, of course, makes 
them relatively useless as density diagnostics, although, 
for He-like ions for which the forbidden line is strong, an 
upper limit can be placed on the density of the X-ray- 
emitting plasma (as well as on the local value of the mean 
intensity). The He-like complexes we observe imply upper 
limits of Ue ~ 10^^ - cm“^. 

It is possible to determine an upper limit for the den¬ 
sity via the Fe XVII I7.I0-to-17.05 A line ratio. The upper 
level of the 17.10 A line is metastable, and so this ratio is 
density sensitive with a critical density of about Ue ~ 10^^ 
cm“^ (Mauche, Liedahl, & Fournier 2001), but the energy 
spacing of this deexciting transition is larger than that for 
the He-like diagnostics, so that photoexcitation requires 
a photon with a wavelength below 400 A. Thus, for all 
practical purposes, in a B star this diagnostic is not sensi¬ 
tive to photospheric UV photoexcitation. In the HETGS 
spectrum of t Sco, we see the 17.10-to-17.05 A ratio in the 
low-density limit. Thus there is no evidence for very dense 
(rie > lO^"* cm“^) X-ray emitting plasma on r Sco. 

5.4. Line Centroids 

In fitting the lines in the grating spectrum, we allowed 
the line centroid wavelength to be a free parameter. These 
results, along with line strengths and widths, are listed in 
Table 1. Although each individual line centroid is roughly 
consistent with the laboratory value of the wavelength, 
when taken as a whole, there is the potential to see an 
aggregate wavelength shift in the ensemble of lines. We 
calculated the weighted mean centroid velocity for the 
roughly two dozen strongest unblended lines and found a 
value of Ucentroid = 14 ± 61 km s“^ (see Fig. 10). The he¬ 
liocentric correction actually shifts this number to the blue 
by about —27 km s“^, making the heliocentric weighted- 
mean line centroid velocity —13 ± 61 km s“^. Thus, there 
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is no direct evidence for redshifted emission that might 
be expected from the shocking of infalling clumps, as in 
the Howk et al. (2000) model, which could be associated 
with the wind component that gives rise to the observed 
redshifted UV absorption. We cannot, however, rule out 
redshifted X-ray emission with net velocities of the order 
of 100 km s“^. 

5.5. Temperature Distribution in the Hot Plasma 

We are deferring detailed DEM modeling of the global 
HETG spectrum to a future paper. The three-temperature 
MEKAL (Mewe, Kaastra, & Liedahl 1995) collisional equi¬ 
librium model that htted the ASCA spectrum of this star 
(Cohen et al. 1997b), does an adequate job of matching 
the gross properties of the HETGS spectrum. Of course, 
with the vast increase in spectral resolution with Chandra, 
this model does not ht in detail. However, it is clear that 
a detailed, global fit to the spectrum will not drastically 
alter the results from the fit to the ASCA spectrum. One 
question left open by the analysis of the medium-resolution 
ASCA data is the extent to which the high-temperature 
component of the spectrum is present, in terms of both 
the overall emission measure of the hot component and 
the actual maximum plasma temperature. 

We find no obvious evidence for plasma with tempera¬ 
tures above the 27 MK lower limit to the hottest tempera¬ 
ture component that was found in the ASCA spectrum, 
and in fact, our detection of the S XV feature, which 
is relatively prominent in the ASCA spectrum, is quite 
marginal in the HETGS spectrum. This is perhaps not 
surprising, however, given the much lower detector effec¬ 
tive area available with HETGS -I- ACIS-S compared to 
ASCA Solid-State Imaging Spectrometer. 

We have made a rough assessment of the temperature 
distribution implied by the H-like to He-like line ratios 
from O, Ne, Mg, and Si. We modeled these ratios using the 
same NLTERT non-LTE ionization/excitation code (Mac- 
Farlane et al. 1993) that we used to model the He-like f/i 
ratios. The best-fit temperature for each of these ions can 
be determined from the observed line ratios under the as¬ 
sumption of coronal equilibrium. We found temperatures 
ranging from less than 3 MK for oxygen to slightly more 
than 10 MK for the silicon. This is in agreement with 
the expectations from the off-the-shelf collisional equilib¬ 
rium codes, such as APEG (Smith et al. 2001), and it 
represents a wider temperature range than is seen in the 
O stars C Pup and C Ori based on this same diagnostic. 
Of course, the presence of a wide range of plasma tem¬ 
peratures is not surprising, and by itself is not a strong 
discriminant among the various physical models of X-ray 
production in hot stars. 

Finally, the wavelength region between 10 A and 12 
A contains several lines from high ionization stages of 
iron, including Fe XXIII at 11.02 A and 11.74 A and 
and Fe XXIV at 10.62 A. The equilibrium temperatures of 
peak emissivity of these lines are roughly 20 MK (APED; 
Smith et al. 2001). These lines are quite strong in the 
spectrum of r Sco. Of the other hot stars observed thus 
far with Chandra HETGS, only 9^ Ori C shows these fea¬ 
tures as strongly as does r Sco. Thus, there is evidence in 
the Chandra spectrum of r Sco for a not-inconsequential 
amount of very hot plasma, qualitatively confirming the 


earlier ASCA results, and demonstrating a similarity be¬ 
tween r Sco and 9^ Ori G that sets both of these young 
hot stars apart from the more standard O supergiants. 

5.6. Time Variability 

Active coronal sources, including pre-main-sequence 
stars, show frequent and strong flaring activity, often in 
conjunction with very high plasma temperatures (e.g., 
AB Dor: Linsky & Gagne 2001; pre-main-sequence stars: 
Montmerle et al. 2002). Wind-shock sources show very 
little variability, and whatever variability is seen has been 
interpreted as being stochastic and reflecting a Poisson 
statistical description of the number of individual sites of 
X-ray emission within the wind (Oskinova et al. 2001). 
Variability levels of 10%, as are seen in the Chandra spec¬ 
tra of T Sco, could then be interpreted as an indication 
that there are roughly 100 individual sites of X-ray emis¬ 
sion at any given time on this star. Of course, arguments 
about characteristic timescales (of, e.g., cooling), as well 
as the evolution of the physical properties of individual 
emission sites, have a bearing on the details of this inter¬ 
pretation. 

Our detection of variability in the oxygen Lya line at 
18.97 A, but in no other line, could indicate that the rela¬ 
tively cooler gas that gives rise to this feature has a more 
variable emission measure than do other lines, and thus 
O VIII might be present in fewer discrete sites within the 
X-ray emitting region, whether it is wind shocks or mag¬ 
netically confined coronal plasma. It could, on the other 
hand, be an indication that the level of absorption by the 
overlying cool wind is variable. The long-wavelength lines 
like O VIII are the most likely lines to be subject to wind 
absorption. However, the expectation is that there is very 
little continuum absorption of the X-rays in a hot star like 
r Sco, which has a relatively low-density wind. We note 
that there is enough signal in several other lines (Si XIII, 
Mg XI, several Fe XVII lines, for example) that variability 
at the level seen in the O VIII line would be detected if 
it were present. We also note that the time dependence 
of the variability in this line tracks the overall X-ray vari¬ 
ability in the spectrum as a whole and also that the even 
longer wavelength lines of O VII and N VII do not have 
enough signal for any variability to be detected. 

6. DISCUSSION 

The qualitative impression one gets from comparing the 
HETGS spectrum of r Sco to that of other stars is that this 
early B star has more in common with late-type coronal 
X-ray sources than with early-type wind sources. This im¬ 
pression is borne out primarily in the analysis of the line 
widths. The emission lines are indeed narrow, although 
they are not as narrow as the same lines seen in Gapella 
and AB Dor. The bulk flow velocities derived from these 
lines, while supersonic, are only 10% - 20% of the wind 
terminal velocity. 

If embedded in a standard CAK wind, the X-ray- 
emitting regions would have to be within just a few tenths 
of a stellar radius from the photosphere based on the ob¬ 
served line widths. However, the He-like f/i ratios put 
constraints on the proximity of the X-ray emitting plasma 
to the photosphere. These results are consistent with a 
height of li?H.-2i?H< (27?* < r < 3i?*), and for Mg XI the 
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plasma must be at least 1.6 i?* above the photosphere. 
There is thus a contradiction, that would seem to imply 
the presence of hot plasma out to a few stellar radii above 
the surface of r Sco, but moving with very small (but 
nonzero) line-of-sight velocity. This picture is inconsis¬ 
tent with standard wind-shock scenarios (which seem to 
apply to O supergiants like C Pup) in which hot plasma is 
embedded within the massive line-driven wind. Our data 
show lines far too narrow for this to be the case. A similar 
situation is seen in the 09.5 star <5 Ori, though the emis¬ 
sion lines in that star are significantly broader than what 
is seen in the spectrum of r Sco (Miller et al. 2002). 

It is also difficult, however, to reconcile these data with 
a generic, solar-type picture of magnetic X-ray activity, 
in which small magnetic loops confine plasma heated by 
magnetic reconnection and the deposition of mechanical 
energy. In this picture, the plasma would be closer to the 
UV-bright photosphere of r Sco than can be accounted for 
by the observed f /i value for magnesium. Furthermore, 
we see no evidence for coronal X-ray flaring in our data. It 
is conceivable that relatively stable, large scale magnetic 
loops could exist on the surface of r Sco, with hot plasma 
confined at the tops of the loops, several stellar radii from 
the surface, but the standard picture of energetic, rela¬ 
tively small spatial scale magnetic activity, as seen in the 
sun, is not consistent with these observations of t Sco. 

One scenario that could explain the two primary results 
from the analysis of the Chandra HETGS spectrum of 
r Sco is the “clump infall” model (Howk et al. 2000), 
which was put forward to explain the hard ASCA spec¬ 
trum and the redshifted UV absorption features seen in 
this star. In this model it is assumed that density en¬ 
hancements (“clumps”) form in the wind of a hot star. 
The authors explored the dynamics of such clumps under 
the combined influence of radiation driving, gravity, and 
wind drag. They found that for the parameters appropri¬ 
ate to r Sco, such density enhancements stall at heights of 
several stellar radii, and fall back toward the star. A bow 
shock would develop around each clump, leading to strong 
shock heating and presumably hard X-ray emission. The 
infalling clumps then also explain the observed redshifted 
UV absorption features. Since the outflowing wind, when 
it passes through the shock front associated with one of 
these density enhancements, would be slowed by a factor 
of 4, this scenario provides a means for explaining the lack 
of strong line-broadening seen in our Chandra spectrum, 
and the heights at which these clumps stall are consistent 
with the He-like f/i ratios. It should be stressed, however, 
that Howk et al. (2000) present no specific mechanism for 
the condensation of these density enhancements out of the 
wind. 

The relative youth of r Sco (Kilian 1994) has suggested 
that a fossil magnetic field may be present in this star. 
The very low projected rotational velocity is also consis¬ 
tent with magnetic spin-down (or, of course, with a pole-on 
viewing angle). Recent modeling of radiation-driven winds 
in the presence of a dipole magnetic field (Babel & Mont- 
merle 1997; Donati et al. 2002; ud-Doula & Owocki 2002) 
has suggested another scenario that may also be relevant 
for T Sco and could explain the formation of density en¬ 
hancements at several stellar radii in this star’s wind. The 
magnetically confined wind shock (MOWS) model predicts 


the confinement of the stellar wind by a strong enough 
dipole field and the subsequent shock heating of wind ma¬ 
terial trapped in the field and forced to collide at the mag¬ 
netic equator. In the case of r Sco, a field of less than 
100 G is required to provide the necessary confinement 
(ud-Doula & Owocki 2002). This model has been success¬ 
fully applied to 0^ Ori C (Donati et al. 2002), which is a 
very young O star with X-ray properties similar to those 
of r Sco. 

Dynamic simulations have shown that plasma in the 
post-shock region in the MOWS model can eventually fall 
back toward the star when it can no longer be supported 
by radiation pressure and magnetic tension (ud-Doula & 
Owocki 2002). By coupling this mechanism to the Howk 
et al. (2000) scenario, we may have a self-consistent, phys¬ 
ical picture of the strong shock heating of wind plasma, 
in a nearly stationary equatorial zone, at a height of sev¬ 
eral stellar radii, with periodic infalling of material. We 
await dynamical simulations of this particular star in order 
to make quantitative comparisons with our observational 
data. In addition, synthesis of X-ray line profiles from the 
dynamical simulations would be very important, since it 
is not clear if the specific observed small line widths can 
be produced within this context. 

The other properties we derived from the r Sco HETGS 
spectrum can perhaps also be explained in this scenario of 
magnetic wind shock heating and confinement combined 
with the infall of density enhancements in the wind. The 
high plasma temperatures are naturally explained by the 
large shock velocities inherent in the collision between the 
fast stellar wind and the equatorial shocked disk (as well 
as with the infalling clumps). 

The lack of very dense plasma argues against very dense 
coronal-type magnetic structures near the photosphere, al¬ 
though some level and type of coronal heating mechanism 
is certainly not ruled out. In fact, recent work (Gharbon- 
neau & MacGregor 2001; MacGregor & Gassinelli 2002) 
suggests that a dynamo mechanism could be operating in 
the interiors of normal OB stars and that if a means of 
bringing the dynamo-generated magnetic fields to the sur¬ 
face could be identified, then solar-type coronal activity on 
hot stars may be quite plausible. The modest time vari¬ 
ability seen in the Chandra spectrum is more significant 
than that seen in O stars that are supposed to be wind 
X-ray sources. This variability level (rms variations at the 
10% level) may simply be indicative of a relatively small 
number of sites of X-ray production on r Sco as compared 
with O stars showing a more standard mode of wind-shock 
X-ray production (Oskinova et al. 2001), and no large flare 
events were seen in this observation of r Sco (or in any 
of the earlier X-ray observations of this star). There¬ 
fore, any coronal model applied to r Sco would need to 
explain the high temperatures (presumably with a strong 
dynamo mechanism of some sort) but also the absence of 
frequent strong flaring. In addition, such a coronal mech¬ 
anism would need to explain the presence of bulk plasma 
velocities of several hundred km s“^, along with at least 
some X-ray-emitting plasma a significant height above the 
surface of the star. 

In summary, the unusual X-ray properties of r Sco, elu¬ 
cidated in great detail for the first time by the superior 
spectral resolution of the Chandra HETGS, now require 
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detailed dynamical modeling in order to test the viability 
of the MOWS model, deduce the role of clump infall, and 
test the alternate scenario of dynamo activity. However, 
we can say with certainty that the standard picture of nu¬ 
merous shocks embedded in a fast stellar wind cannot ex¬ 
plain the narrow X-ray lines and the formation heights of 
several stellar radii of the very hot X-ray emitting plasma 
on this young hot star. 


We wish to thank Eric Levy, Carolin Cardamone, James 
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1089A to Swarthmore Gollege, Prism Gomputational Sci¬ 
ences, and the University of Wisconsin. Work was per¬ 
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ergy by the University of Galifornia Lawrence Livermore 
National Laboratory under contract No. W-7405-Eng-48. 
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Fig. 3.— X-ray light curve formed from the combined MEG H-l and -1 order counts, with 1000 second bins. The mean count rate is 
indicated by the dashed line. The hypothesis of a constant source can be rejected at a more than 99.99% confidence level. 
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Fig. 4.— X-ray light curve formed from the combined MEG -hi and -1 order counts in the O VIII line at 18.97 A, with 10,000 s bins. The 
mean count rate is indicated by the dashed line. The hypothesis of a constant source can be rejected at a more than 99% confidence level. 
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Fig. 5.— MEG +1 and -1 order observation of the neon Lya line (histogram) with an intrinsically narrow model (convolved with the 
instrument response). The fit to the data show a statistically significant line width. 



Fig. 6.— Derived line widths (HWHM) for three strong lines in seven stars: two stars representative of coronal sources (Capella and AB 
Dor: open symbols connected by dotted lines), r Sco (filled diamonds and solid line), and four O stars (filled symbols and dashed lines), which 
are presumably wind X-ray sources. 
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Fig. 7.— HWHM for the strongest unblended lines in the HETGS spectrum of r Sco. The best-fit linear function is indicated by the solid 
line with the 95 % confidence limits shown as dashed lines. 
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Fig. 8.— MEG H-l and -1 order data for the helium complexes of S XV, Si XIII, Mg XI, Ne IX, and O VII along with the best-fit models. 
Note that the error bars shown here (Sherpa’s “chi Gehrels” option) represent Poisson noise plus a constant. We used the Cash C statistic as 
our goodness-of-fit indicator, and this statistic does not make use of the formal errors. 
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Fig. 9.— Calculation of the density and mean-intensity sensitivity of the f/i line ratio for Mg XL Note that the sensitivity to the mean 
intensity enters via the distance from the photosphere. The distance of the plasma from the photosphere is indicated for each model in units 
of the stellar radius. The range from the data is indicated by the shaded area. 
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Fig. 10.— Line centroids derived from fitting Gaussian models to the strongest lines in the HETGS spectrum. The dashed line indicates 
the weighted mean centroid velocity. 
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Table 1 


Emission Lines in the HETGS Spectrum 


Ion 

Mab 

-^obs 

Flux 

Half-Width“ 

Centroid Velocity^ 

log Tmax —emiss 


(A) 

(A) 

(10“^ ph s“^ cm“^) 

(km s“^) 

(km s“^) 

(K) 

sxv-* 

5.0387 ± 0.002 


0.936+;^™ 



7.1 

SKV-* 

5.0631 ± 0.002 


0.464+;?!^ 



7.1 

SXV* 

5.1015 ± 0.002 


0.944+;^®« 



7.1 

Si XIV 

6.1821 ± 0.0003 

6.1832 ± 0.0011 


4241?! 

c-i+53+15 

^^-53-15 

7.2 

Si XIII 

6.6479 ± 0.002 

6.6495 ± 0.0006 

3.771 

2751“ 

70 + 27+90 
' ^-27-90 

7.0 

Si XIII 

6.6882 ± 0.002 

6.6881 ± 0.0014 

lA39t\tl 

4221™® 

4+63+90 

‘*-63-90 

7.0 

Si XIII 

6.7403 ± 0.002 

6.7376 ± 0.001 

3.0l4t;^^® 

3971“ 

-I on+^5 + 89 
-*-■^^-45-89 

7.0 

Fe XVIII 

7.0889 ± n/a 

7.0952 ± 0.006 

0.378t;J^? 

169+211 

10^-169 

266l??"‘ 

7.2 

Mg XII 

7.1062 ± 0.00004 

7 11 

' • ' —0.0032 


^au_i90 

1 qo+®5+2 
*^^-135-2 

7.0 

Mg XI 

7.8503 ± 0.002 

7.8477 ± 0.0016 

0.827j:;lJ'; 

q-1 -^+106 

_99+61 + 76 
^^-61-76 

6.8 

Al XII 

7.8721 ± 0.003 

7.8814 ± 0.0032 

0.339j:;o*4 

10=^-189 

^^4+114+122 

'^'J^-114-122 

6.9 

Mg XII 

8.4210 ± 0.0001 

8.422 ± 0.001 

2.802tit°i 

3031“ 

07+36+2 
' -36-2 

7.0 

Mg XI 

9.1688 ± 0.003 

9.173 ± 0.001 

A.230tfjl 

3061™ 

-I 07+29+98 

AO < _29-98 

6.8 

Mg XI 

9.2308 ± 0.003 

9.2334 ± 0.0012 

2.634tlll 

3171“ 

00+39 + 98 

00-39-98 

6.8 

Mg XI 

9.3143 ± 0.003 

9.3124 ± 0.002 

1.263-t-\ll 

4281?? 

0-1+64+97 

'JA_g4_97 

6.8 

FeXXI 

9.4797 ± n/a 

9.4762 ± 0.0014 

i.soetlfo 

2631?! 

_lll+44 

A-LA_44 

7.0 

Ne X 

10.2388 ± 0.00005 

10.239 ± 0.002 

2.069t;^39 

2941?! 

c: + 59+l 
^-59-1 

6.8 

Fe XVIII 

10.3603 ± n/a 

10.359 ± 0.002 

1.15lj::^t4 

Q-i -1 +123 
■^-l-'--171 

-94l?? 

6.9 

Fe XVIII 

10.537 ± n/a 

10.536;“;™"2 

0.8151:^34 

39l!™ 

-29l?! 

6.9 

Fe XXIV 

10.619 ± 0.02 

10.625 ± 0.004 

2.057l;|3^ 

482+151 

.,70 + 113+565 

A'U-113-565 

7.3 

Fe XVII 

10.77 ± 0.002 

10.773 ± 0.002 

1.850l;|t3 

-‘-^■^-108 

0 .+57+56 
^4-57-56 

6.8 

Fe XVII 

11.254 ± 0.002 

11.251 ± 0.001 

2.754till 

2731“ 

00+27+53 

^^-27-53 

6.8 

Fe XVIII 

11.326 ± 0.004 

11.321 ± 0.002 

2.975l;|™ 

4671?? 

qo+53 + 106 
AO^-53-106 

6.9 

Ne IX 

11.544 ± 0.01 

11.546 ± 0.004 

3.18lj:;339 

5091?™ 

ro + 104+260 
^■^-104-260 

6.6 

Fe XXIII 

11.736 ± 0.004 

11.745 ± 0.002 

2.939tllt 

000 +113 

oqn+51 + 102 
■^OU-5i_io2 

7.2 

Fe XXII 

11.770 ± 0.003 

1 1 77c:-l-u.uUJ 

11. 1 <0_Q QQ2 

l.696l;«^ 

67l!!® 

107+76 + 76 
' -51-76 

7.1 

Ne X 

12.134 ± 0.0001 

12.132 ± 0.001 

16.8321'™® 

2001™ 

_46+25+i 

^^^-25-1 

6.8 

FeXXI 

12.284 ± 0.002 

12.275 ± 0.001 

9A14tlll 

4331?? 

000+49 + 24 
■^■iU-49-24 

7.0 

FeXXI 

12.393 ± 0.005 

12.401 ± 0.004 


198l??t 

194+97+121 

a»^_97_121 

7.0 

FeXX 

12.823 ± 0.005 

12.833 ± 0.002 

7.847l;®“ 

5051?! 

0^0+47+117 

7.0 

Ne IX 

13.447 ± 0.004 

13.450 ± 0.001 

10.085l'|™ 

1781?? 

00+22 + 89 
'^^-22-89 

6.6 

FeXIX 

13.518 ± 0.002 

13.514 ± 0.002 

8.2611;™® 

312I?! 

00+44+44 

oy_44_44 

6.9 

Ne IX 

13.553 ± 0.005 

13.553 ± 0.001 

10.850l'J®® 

1 

-*-^'-35 

0 + 22+111 
'-’-22-111 

6.6 

FeXIX 

13.795 ± 0.005 

13.782 ± 0.002 

2.394l;®“| 

1141??? 

Ooq+44+109 

■^00_44_ j^Q9 

6.9 

Fe XVII 

13.825 ± 0.002 

13.824;“;™^ 

3.4961;®™ 

3281?™ 

00+87+43 

^^-65-43 

6.8 

Fe XVIII 

14.208 ± 0.003 

14.203 ± 0.001 

10.4551 

2841?? 

_106+21 + 63 
J-UU-21-63 

6.9 

Fe XVIII 

14.256 ± 0.005 

14.255 ± 0.002 

4.68i1;®33 

3141?? 

01+42 + 105 
■^-'--42-105 

6.9 

Fe XVIII 

14.534 ± 0.003 

14.536 ± 0.002 

5.4851;®™ 

2041?? 

.1+41 + 62 
^-*--41-62 

6.9 

Fe XVIII 

14.571 ± 0.011 

14.575 ± 0.004 

l'215l;«® 

00 + 233 
^^-29 

00+82 + 226 
0-^-82-226 

6.9 

Fe XVII 

15.014 ± 0.001 

15.017 ± 0.001 

32.36211®™ 

2281?? 

00+20 + 20 
'-’'-’-20-20 

6.7 

Fe XVII 

15.261 ± 0.002 

15.264 ± 0.001 

13.100+1“™ 

226l|? 

+20 + 39 
'A =^-20-39 

6.7 

Fe XVIII 

15.625 ± 0.003 

15.628 ± 0.002 

5.1511;™? 

1761?? 

c-o+38 + 58 
'JO-38-58 

6.7 

O VIII 

16.006 ± 0.0001 

16.006 ± 0.001 

13.47911,““! 

1491?! 

1 q+19 + 37 
-''^-19-37 

6.5 

Fe XVIII 

16.071 ± 0.003 

16.075 ± 0.001 

8.0351;?™ 

1291?? 

7^ + 19 + 56 
'’*-19-56 

6.8 

Fe XVII 

16.780 ± 0.002 

16.778 ± 0.001 

fjq q'7-1+1-358 

ZO.OI -L_i.309 

1481?? 

qo + 18 + 36 
'^'*-18-36 

6.7 

Fe XVII 

17.051 ± 0.001 

17.054 ± 0.001 

29.20711,“® 

2131?? 

c.q+18 + 18 
^'*-18-18 

6.7 

Fe XVII 

17.096 ± 0.001 

17.099 ± 0.001 

2fi 1 14+1-534 

130l?? 

i-q+18 + 18 
^'*-18-18 

6.7 

O VIII 

18.969 ± 0.0001 

18.970 ± 0.001 

50.569111™ 

29ll?! 

10 + 16 + 1 
-*-®-16-l 

6.5 

O VII 

21.602 ± 0.007 

21.604 ± 0.003 

01 710+2-708 

Zl. 1 -1-^_2.506 

300l?! 

q-+42 + 97 
"*^-42-97 

6.3 

O VII 

21.804 ± 0.007 

21.809 ± 0.002 

07 QY1+5-110 
z / .z / -•-_2 g9Q 

177+55 
‘ ‘ -32 

74+28 + 96 
'‘*-28-96 

6.3 

N VII 

24.781 ± 0.0001 

24.783 ± 0.003 

22.01911??? 

3091?? 

Or+36+1 

^^-36-1 

6.3 


“Derived HWHM. 

^The first quoted uncertainties are statistical and based on the fit of the Gaussian model, while the second quoted uncertainties are 
systematic and are based on the uncertainties in the laboratory rest wavelengths as compiled in APED. 

‘^Temperature of maximum emissivity; from APED. 

‘^Because the position and width of the Gaussian models were fixed for these lines, we do not list best-fit values or uncertainties for 
these quantities. 
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Table 2 

Emission-Line Blends 


Ion 

\ a 

Mab 

Flux 


(A) 

(10“^ photons s“^ cni“^) 

FeXX 

10.00 

q f^qo-t-U.:^»D 
o.uoo_Q 263 

Ni XIX 

10.58 

1 SQQ+0-386 

-i-.oyy_o.342 

Fe XVIII, Fe XIX, Fe XXIV 

10.65 

U.405_o.i7g 

Fe XXIII 

11.00 

0 qi^7+0-384 

Z.OU/ _0.35i 

FeXXI 

12.42 

A rrr+0.463 
^.000_o.436 


“Approximate mean wavelength of the blended feature. 


Table 3 

Helium-like Forbidden-to-Intercombination Ratios 


Ion 


f/i 

Formation Radius 


(A) 

{R*) 

0 VII 

1623 

< 0.13 

< 10 

Ne IX 

1263 

< 0.03 

< 2 

Mg XI 

1036 

0.48 ±0.07 

2.6-2.9 

Si XIII 

865 

2.101-28 

1.1 - 1.5 

SXV 

743 

2.03418 

< 20 


Wavelength of the ® P transition that depletes the upper level of the forbidden line and enhances the intercombination line. 



